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ABSTRACT 

We undertook a time-series photometric multi-site campaign for the rapidly oscil- 
lating Ap star HD 99563 and also acquired mean light observations over two seasons. 
The pulsations of the star, that show flatter light maxima than minima, can be de- 
scribed with a frequency quintuplet centred on 1557.653 /iHz and some first harmonics 
of these. The amplitude of the pulsation is modulated with the rotation period of the 
star that we determine with 2.91179 ± 0.00007 d from the analysis of the stellar pul- 
sation spectrum and of the mean light data. We break the distorted oscillation mode 
up into its pure spherical harmonic components and find it is dominated by the I = 1 
pulsation, and also has a notable t — 3 contribution, with weak t — and 2 com- 
ponents. The geometrical configuration of the star allows one to see both pulsation 
poles for about the same amount of time; HD 99563 is only the fourth roAp star for 
which both pulsation poles are seen and only the third where the distortion of the 
pulsation modes was modelled. We point out that HD 99563 is very similar to the 
well-studied roAp star HR 3831. Finally, we note that the visual companion of HD 
99563 is located in the S Scuti instability strip and may thus show pulsation. We show 
that if the companion was physical, the roAp star would be a 2.03 M Q object, seen 
at a rotational inclination of 44°, which then predicts a magnetic obliquity (3 = 86.4°. 

Key words: stars: variables: other - stars: oscillations - techniques: photometric - 
stars: individual: HD 99563 - stars: individual: XY Crt 



1 INTRODUCTION 

Where the lower part of the classical instability strip inter- 
sects the main sequence, three distinct classes of multiperi- 
odically pulsating variables can be found. The 7 Doradus 
stars pulsate in gravity modes of high radial order and have 
periods of the order of one day (Kaye et al. 1999). The S 
Scuti stars have periods of the order of a few hours (Breger 
1979) and are thus pressure and gravity mode pulsators of 
low radial order. 

The fastest pulsations in this domain in the HR diagram 



are however excited in the rapidly oscillating Ap (roAp) 
stars (Kurtz 1982, Kurtz & Martinez 2000), with typical 
periods around 10 minutes, indicating pressure modes of 
high radial overtones. The photometric semi-amplitudes as- 
sociated with these pulsations are in most cases only a few 
mmag, which makes them difficult to detect. The roAp stars 
are also remarkable due to their spectral peculiarities, since 
they are pulsating representatives of cool magnetic Ap stars 
of the SrCrEu subtype. 

The magnetic fields of Ap stars cause elemental segre- 
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gation on the stellar surface. In other words, the chemical 
elements are arranged in patches on the stellar surface (see, 
e.g., Kochukhov et al. 2004 or Lueftinger et al. 2003 for well- 
documented examples). This causes modulation of the mean 
apparent brightness of the star with the rotation period. As 
the chemical elements show some alignment with the mag- 
netic poles, the rotational light curves show a single-wave 
structure if only one magnetic pole is seen and a double- 
wave variation if both magnetic poles come into view during 
a rotation cycle. 

The strong magnetic fields present in the roAp stars 
are usually not aligned with the stellar rotation axis. Since 
the pulsation axis does however coincide with the magnetic 
axis, the roAp stars are oblique pulsators (Kurtz 1982). This 
means that the pulsation modes excited in the roAp stars 
are seen at different aspects during the rotation cycle, which, 
for nonradial modes, causes amplitude modulation over the 
rotation period (which can then be inferred) and allows us 
to put constraints on the geometry of the pulsator. 

The oblique pulsator model (OPM), in its simplest 
form, then predicts that (except for special geometric ori- 
entations of the axes) dipole (£ — 1) pulsation modes will be 
split into equally spaced triplets. Quadrupolc (£ = 2) modes 
will give rise to equally spaced frequency quintuplets, where 
the spacing of consecutive multiplet components is exactly 
the stellar rotational frequency. 

The magnetic field of the roAp stars affects the pul- 
sations in two additional ways. Firstly, the pulsation fre- 
quencies are shifted with respect to their unperturbed value 
(Cunha & Gough 2000) and secondly, the pulsation modes 
are distorted so that a single spherical harmonic can no 
longer describe them fully (e.g., Kurtz, Kanaan & Martinez 
1993, Takata & Shibahashi 1995). This is observationally 
manifested by the presence of additional multiplet compo- 
nents surrounding the first-order singlets, triplets or quintu- 
plets, which are spaced by integer multiples of the rotation 
frequency. 

For a long time, the predominant observing method 
for studying the pulsations of roAp stars was time-resolved 
high-speed photometry, yielding interesting results on the 
pulsational behaviour, geometry and asteroseismology of 
these stars (e.g., see Matthews, Kurtz & Martinez 1999). 
The most recent observational advances however came from 
time-resolved spectroscopy: since the vertical wavelengths 
of the pulsation modes are comparable to the size of the 
line-forming regions in the atmospheres of these stars, 
the vertical structure of the atmospheres can be resolved 
(Ryabchikova et al. 2002, Kurtz, Elkin & Mathys 2003). 

Spectroscopy is also more sensitive to the detection of 
roAp star pulsations than photometry. Due to the vertical 
stratification of chemical elements in their atmospheres (e.g. 
Ryabchikova et al. 2002) the pulsational radial velocity am- 
plitudes of some spectral lines (most notably rare earth ele- 
ments) can reach several kilometres per second. 

The most extreme example of such high radial velocity 
amplitudes is HD 99563 (XY Crt). This star was photometri- 
cally discovered to pulsate by Dorokhova & Dorokhov (1998) 
and was confirmed by Handler & Paunzen (1999). Elkin, 
Kurtz & Mathys (2005) discovered pulsational radial veloc- 
ity variations with semi-amplitudes of up 5 kins -1 for some 
EuII and Tmll lines in their time-resolved spectroscopy of 
this star. Such high amplitudes are capable of yielding in- 



formation on the structure of the atmospheres of Ap stars 
with the best possible signal to noise. 

However, one important piece of information that spec- 
troscopy cannot supply at this point is detailed knowledge of 
the stellar pulsation spectrum. The reason is that the largest 
telescopes are necessary for obtaining spectra of the required 
time resolution and signal to noise. However, observing time 
on these telescopes is sparse. Therefore, lengthy photomet- 
ric measurements of these stars on small telescopes are still 
required to decipher the pulsational spectra fully. 

To this end, we included HD 99563 as a secondary target 
in a multi-site campaign originally devoted to the roAp star 
HD 122970 (Handler et al. 2002) to be observed at times 
when the latter star was not yet accessible. However, HD 
99563 turned out to be quite interesting, which is why we 
continued to observe it after the original campaign was fin- 
ished. We also acquired some mean-light observations of HD 
99563 in an attempt to determine its rotation period. In this 
paper, we report the results of these measurements. 



2 OBSERVATIONS AND REDUCTIONS 
2.1 The high-speed photometry 

Our multi-site time-series photometric observations were ob- 
tained with seven different telescopes at four observatories: 
the 0.75m T6 Automatic Photometric Telescope at Fairborn 
Observatory (FAPT) in Arizona, the 0.5m, 0.75m, 1.0m and 
1.9m telescopes at the South African Astronomical Observa- 
tory (SAAO), the 0.6m reflector at Siding Spring Observa- 
tory (SSO) in Australia and the 0.9m telescope at Observa- 
torio de Sierra Nevada (OSN) in Spain. Whereas the latter 
observations consisted of Stromgren uvby photometry, the 
other telescopes used the B filter only. The measurements 
are summarised in Table 1. 

The time-series photometry was generally acquired as 
continuous 10-second integrations. Large apertures (> 35 
arcseconds) were employed to minimise the contributions of 
seeing and guiding. The observations of the target star were 
interrupted at irregular intervals (depending on the bright- 
ness and position of the Moon) to obtain measurements of 
sky background. 

Reductions were performed by first correcting for coinci- 
dence losses (dead time correction), then for sky background 
and extinction. Some low-frequency filtering to remove resid- 
ual transparency variations or tube sensitivity drifts with 
time scales longer than 30 minutes was applied by fitting 
low-order polynomials to the data and by subtracting these. 
Finally the measurements were binned into 40-second in- 
tegrations, the times of the observation were converted to 
Heliocentric Julian Date (HJD) and all the data were joined 
into a combined single light curve. 

Some of our nightly measurements are shown in Fig. 1, 
clearly containing some interesting features. Firstly, the pul- 
sation amplitude can vary rather dramatically from almost 
zero to more than 10 mmag peak-to-peak, and it can also 
be seen to vary during a single night. Secondly, the shape of 
the light curves near maximum amplitude is peculiar, with 
the light maxima being flatter than the light minima. 
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Figure 1. Some light curves of HD 99563 from our multi-site observations showing its diverse nightly behaviour. The solid lines represent 
the multifrequency fit to be listed in Table 2, to guide the eye. 



2.2 The mean light observations 

Differential multicolour photometry of HD 99563 was ob- 
tained with the 0.75m Automatic Photometric Telescope at 
SAAO (Martinez et al. 2002) from March 2001 until June 
2004. Single-channel differential measurements with respect 
to the single comparison star HD 99506, which is not known 
to be variable, in the Johnson-Cousins UBVRI filters were 
acquired. Again, large apertures (> 35") were used to elim- 
inate effects of scintillation and seeing. 

Two nightly UBVRI measurements of HD 99563 were 
usually taken in between two measurements of HD 99506, 
ensuring that variations in sky transparency can be readily 
compensated for, but also aiming at averaging out the effects 
of the rapid oscillations on the mean magnitude. As the two 
integrations on the target were separated by about 380 s, 



about half the oscillation period, the pulsational signal was 
strongly suppressed. 

The data of both stars were reduced by correction for 
coincidence losses, sky background and extinction, and the 
timings were heliocentrically corrected. Differential magni- 
tudes between the variable and the comparison star were 
computed by simple linear interpolation. After rejection of 
outliers and merging of the two nightly measurements, a to- 
tal of 150 differential UBVRI magnitudes of HD 99563 were 
available. 



3 FREQUENCY ANALYSIS (I) 
3.1 The high-speed photometry 

The reduced B filter light curves were analysed with 
the computer programme Period98 using single-frequency 
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Table 1. Journal of the observations of HD 99563; JV40 is 
the number of 40-second data bins obtained in the respective 
night and <r is the rms scatter of the residual light curves after 
prewhitening and low-frequency filtering; the remaining columns 
are self-explanatory. 



Telescope Date Start Length A^4q a 

(UT) (UT) (hr) (mmag) 



Total B 125.49 9728 1.82 

uvby 6.05 444 7.7/3.9/3.7/4.3 



Fourier and multiple-frequency least-squares techniques 
(Sperl 1998). This programme has, amongst others, the ca- 
pabilities of calculating Fourier periodograms, calculating 
the best fit to the data with all given frequencies simultane- 
ously by minimising the residuals between light curve and 
fit. Within the multifrequency fitting, combination frequen- 
cies (that are in our case expected due to the shape of the 
light curves) can be fixed to their expected values and their 
amplitudes and phases can be optimised together with the 
other independent parameters. 

We show the amplitude spectrum of the combined light 
curve in the uppermost panel of Fig. 2. Two main regions 
of power are visible, one centered around 1500 /iHz and a 
weaker one around 3100 /iHz. The first of the sequence of 
panels on the left-hand side of Fig. 2 shows the spectral 
window of the data at an extended scale, computed as the 
Fourier transform of a single noise free sinusoid with a fre- 
quency of 1553.67 /iHz and an amplitude of 2.7 mmag. We 
see some alias patterns about 34.7 /iHz away from the cen- 
tral peak. This is the 3 cd~ alias, which is most prominent 
in our data because of the longitude distribution of the ob- 
servatories used combined with the short duration of most 
of the runs. Fortunately, neither this alias, nor any other, 
will have any influence on our results. 

The second panel on the left-hand side contains the am- 
plitude spectrum of the data. It is clearly more complicated 
than the spectral window. Prewhitening the strongest sig- 
nal from the data we end up with the residual amplitude 
spectrum in the third (middle) panel of Fig. 2, left-hand 
side. Further prewhitening of the strongest signal there re- 
veals another one, and continuing this procedure results in 
the detection of a fourth frequency in this domain, where no 
more signals can now be detected. 

We continue the frequency analysis in the region around 
3100 /iHz, showing it in the sequence of panels on the right- 
hand side of Fig. 2. Following the same strategy as before, we 
detect five signals in this frequency domain. After prewhiten- 
ing these, we arrive at the residual amplitude spectrum in 
the lowest panel of Fig. 2. The noise level at frequencies 
below ~ 800/iHz is artificially decreased due to our low- 
frequency filtering. 

Consequently, we have separately analysed our longest 
(AT > 2 hr) light curves from the best nights, leaving them 
unfiltered for low frequencies, and searched them for period- 
icities. We found no such variations, with detection thresh- 
olds of 1.8 mmag for frequencies below 200 /iHz, 0.8 mmag 
for frequencies between 200 and 450 /iHz, and 0.5 mmag for 
frequencies between 450 and 1000 /iHz, respectively. 

More interesting is the excess mound of amplitude at 
frequencies around 1500 /iHz on top of the essentially mono- 
tonically decreasing noise that can be discerned in the low- 
est panel of Fig. 2. This is the region where four signals 
have already been detected. The excess mound of amplitude 
may indicated the presence of further signals or of ampli- 
tude/frequency variability of some of the periodicities de- 
tected previously. In any case, a preliminary frequency so- 
lution derived in this way is listed in Table 2. 

We have only accepted signals that satisfy the statis- 
tical criterion of S/N > 4 for independent oscillations and 
S/N > 3.5 for combination frequencies. These variations are 
regarded as statistically significantly detected. We refer to 
Breger et al. (1999) for a more in-depth discussion of this 
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Figure 2. Amplitude spectra of HD 99563 from our multi-site observations with successive prewhitcning of the detected frequencies. 



criterion, that has turned out to be extremely reliable in the 
past. 

The independent frequencies in Table 2 can, at first 
sight, not be unambiguously be interpreted within the OPM. 
Although /a, /i and fa form a triplet with a spacing of about 
7.95 /xHz, /i , fi and /2 can also be reconciled with a triplet 



with a spacing of approximately 3.98 /iHz, allowing for a 
~ 3cr shift in frequency f^. In this context it is important to 
note that the formal error estimates we use here are believed 
to underestimate the real errors by about a factor of 2 (Han- 
dler et al. 2000, Jerzykiewicz et al. 2005), and therefore we 
cannot simply reject the latter interpretation. Consequently, 
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Table 2. Preliminary multifrequency solution for our B filter data 
of the roAp star HD 99563, with the frequencies around 1500 fiHz 
left as free parameters. Error estimates correspond to Icr values 
are determined following Montgomery & O'Donoghue (1999). 



Table 3. The amplitudes and phases of the mean light variations 
of HD 99563, once more with error estimates from Montgomery 
& O'Donoghue (1999) The phases are with respect to pulsation 
amplitude maximum at HJD 2452031.29627. 



ID 



Frequency 
(fiHz) 



B Ampl. 
(mmag) 



S/N 
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Figure 3. Amplitude spectra of the mean light observations of 
HD 99563 in the U filter. Upper panel: spectral window func- 
tion of the time series. Middle panel: amplitude spectrum of the 
data themselves. Lower panel: residual amplitude spectrum after 
prcwhitening the strongest signal from the data. Note the five 
times larger ordinate scale of the lower panel. 



all four frequencies could even be part of an equally spaced 
quintuplet with a component missing near 1549.7 ^iHz. 



3.2 The mean light observations 

We again used Period98 for the frequency analysis. Ampli- 
tude spectra of the measurements in the five UBVRI niters 
were computed. Although the nominal Nyquist frequency of 
our observations is 0.5 cd _1 (~ 6 /kHz), we carefully exam- 
ined a frequency range up to 25 pHz. We found variability 
in all the five passbands, with those in U having the highest 
amplitude, as shown in Fig. 3. 

The upper panel of this graph contains the spectral 
window function of the data, calculated as the amplitude 
spectrum of a single noise-free sinusoid with a frequency of 



Filter 


Amplitude 


Phase 




(mmag) 


(rad) 


U 


31.3 ± 1.0 


-0.06 ±0.03 


B 


20.4 ± 0.9 


-0.00 ±0.04 


V 


4.6 ± 0.8 


+3.09 ± 0.17 


R 


9.2 ± 0.8 


+3.09 ± 0.09 


I 


7.3 ± 0.7 


+3.11 ± 0.10 



7.95 /jHz and an amplitude of 31 mmag. There is strong 
1 cd -1 aliasing (the peak near 20 /iHz), and also reflection 
of the aliasing pattern at frequency zero (peaks near 4 and 
15 piHz). Nevertheless, the input signal produces the tallest 
peak in the window function. This is why we can be sure 
that the peak labelled in the amplitude spectrum of the data 
(middle panel of Fig. 3) is also the correct frequency present 
in the measurements, which is confirmed by the prewhitened 
amplitude spectrum in the lower panel of this figure that 
contains noise only. 

We encountered the same situation in the B filter data, 
where the mean light variations also have considerable am- 
plitude. In VRI the situation is more complex, as the S/N 
of the variability is lower. Significant peaks at the same fre- 
quencies are however still present, but due to the effects of 
noise the tallest maxima are found at alias frequencies in 
V and R. After prewhitening the signal near 7.95 ^Hz, the 
residual amplitude spectra in all filters contain noise only; 
in particular, no evidence for harmonics or subharmonics of 
this frequency can be detected. 

Consequently, we can describe the variability present 
in the mean light variations of HD 99563 with a single fre- 
quency of 7.9498 ± 0.0002 /iHz, as determined from a non- 
linear least-squares fit to the U data, which have the best 
signal-to-noise. This frequency is consistent with one of the 
possible multiplet spacings within the pulsational signals as 
determined in the previous section, but the ambiguity of 
whether this frequency or half of it corresponds to the rota- 
tion period of HD 99563 is not yet resolved. The amplitudes 
and phases of the mean light variations with respect to the 
frequency of 7.9498 /iHz are listed in Table 3. 

This table shows that the mean light maxima in the 
U and B filters occur at the same time as pulsation am- 
plitude maximum in the light curves. In the V, R and I 
filters, mean light minimum coincides with pulsation am- 
plitude maximum. We find no statistically significant phase 
lag of the mean light extrema relative to pulsation amplitude 
maximum. 



4 THE ROTATION PERIOD OF HD 99563 

To infer the pulsational and magnetic geometry of HD 
99563, its rotation period has to be known with certainty. 
However, at this point it is not clear whether the 7.95 /iHz 
signal in the mean light variations and the frequency split- 
ting within the pulsational signals correspond to the rotation 
period or to half the rotation period. 
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Figure 4. Mean UBVRI light variations of HD 99563 phased 
with a period of 2.91179 d. One and a half rotation cycles are 
shown. The solid line is the single-frequency fit from Table 3. 



Since the magnetic field of the star also varies with as- 
pect, we can invoke published magnetic field measurements 
for HD 99563 to resolve this ambiguity. Hubrig et al. (2004) 
measured a longitudinal magnetic field strength of —688 ± 
145 G on HJD 2452494.479, whereas Kudryavtsev & Ro- 
manyuk (2005) measured +580±100 G on HJD 2453395.550. 

Assuming a rotation frequency of 7.9498 ± 0.0002 fiHz, 
the two measurements would have been taken 618.91 ± 0.02 
rotational periods apart. Since these measurements would 
then have been taken at nearly the same rotation phase, 
but the magnetic field showed a reversal between these two 
measurements, such a rotation frequency can be ruled out. 
Using half this value (3.9749 ± 0.0001 ^Hz) implies that 
309.455 ± 0.009 rotations have gone by, which is perfectly 
consistent with the magnetic field reversal between the two 
magnetic field observations and the absolute magnetic field 
strengths measured. 

The correct value of the rotation period of HD 99563 is 
therefore close to 2.91 d. We show the corresponding phase 
diagram of the mean light variations of the star in Fig. 4. 
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Figure 5. The frequency quintuplet detected in the light curves 
of HD 99563. 



OPM. We first examine the variations in the frequency do- 
main near 1550 ^iHz. The spacings within all the frequencies 
already detected in this range are consistent with this rota- 
tion frequency and thus belong to the same multiplet. We 
therefore searched for possible further multiplet components 
in the residual amplitude spectra that have S/N > 3.5. We 
indeed detected a peak near /s = 1569.6 fjMz that is con- 
sistent with this splitting and has S/N > 3.6. Adding this 
signal to the four previously detected, we find a symmetri- 
cal series of peaks centred on f± = 1557.6 /iHz, for which we 
show a schematic graphical representation in Fig. 5. 

We can now derive a final frequency solution. To this 
end, we have started with the frequencies from the prelimi- 
nary solution in Table 2, added /s and fitted these ten signals 
to the data. During this process, we fixed the signal frequen- 
cies to equal splitting, using a starting value of 3.9749 /iHz, 
and we also demanded that harmonics occurred at exactly 
twice the "parent" frequency within Period98. With this 
procedure, we do not only determine final frequency values, 
but we also reduce the number of free parameters in the fit 
to a minimum. This should guarantee the most reliable and 
most stable frequency solution, and the rotation frequency 
will also be determined most accurately that way. The final 
result of our frequency analysis can be found in Table 4. 

The splitting between the rotationally split components 
of the pulsation mode of HD 99563 is 3.9751 ± 0.0003, which 
is the same as half the frequency of the mean light variations 
within the errors. Since we know that this splitting corre- 
sponds to the rotation frequency, we can fix the latter by 
averaging the two results into 3.9749 ± 0.0001 ^iHz, corre- 
sponding to a rotation period of 2.91179 ± 0.00007 d. 

The referee has noted that relative phases of the first- 
order combination frequencies (0(2^) — 2(j)(v) — 0.01 ± 0.20; 
0(2z/_i) - 20(z/_i) = -0.27 ± 0.10; 4>{2v +1 ) - 2<j>(v +1 ) = 
0.14 ± 0.16) are consistent with zero. This is an indication 
that the frequencies around 3115 /iHz are indeed the first 
harmonics of the main pulsational signals near 1557 pHz. 



5 FREQUENCY ANALYSIS (II) 

Knowing that the rotation frequency of HD 99563 is near 
3.975 pHz, we can now determine a final multifrequency so- 
lution for the pulsation data under the assumption of the 



6 DISCUSSION 

The pulsational light variations of HD 99563 are due to a 
single pulsation mode at a frequency of 1557.6530 /iHz. Due 
to the rotation of the star with a period of 2.91179 d, we 
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Table 4. Final multifrcqucncy solution for our B filter data of the 
roAp star HD 99563, again with error estimates from Montgomery 
& O'Donoghuc (1999; errors on frequencies and amplitudes are 
given in Table 2). Pulsation phases are with respect to pulsation 
amplitude maximum at HJD 2452031.29627. 



ID 






Frequency 
OuHz) 


B Ampl. 

(mmag) 


Phase 
(rad) 


v = 


U 




1557.6530 


0.29 


+2.63 ±0.09 


V-l 


= h- 


frot 


1553.6779 


2.67 


+2.87 ±0.01 


"+1 


= h + 


frot 


1561.6281 


1.82 


±2.83 ±0.01 




= / 4 - 


3 frot 


1545.7276 


0.45 


+2.82 + 0.06 


^+3 


= U + 


3 frot 


1569.5784 


0.18 


±3.00 ±0.15 


2v 






3115.3060 


0.32 


-1.03 + 0.08 


2v- 


i 




3107.3557 


0.25 


-0.27 + 0.10 


2v + 


i 




3123.2562 


0.17 


-0.76 + 0.16 


2v -4 


" frot 




3119.2811 


0.15 


-0.51 + 0.18 


2u - 


- if ro t 




3099.4055 


0.13 


+0.55 + 0.21 


frot 






3.9751 







see it at different aspect during the rotation cycle, which 
gives rise to an equally spaced frequency triplet. Two more 
multiplet components separated by three times the rotation 
frequency were also detected. These are an indication that 
the pulsation mode of HD 99563 is not a pure dipole, but is 
distorted by the effect of its magnetic field. 

To examine the nature of this pulsation mode, we first 
need to know the inclination of the stellar rotation axis to 
the line of sight, which in turn requires an estimate of the 
stellar radius. Consequently, we first infer the star's position 
in the HR Diagram. 

HD 99563 is a visual binary with a secondary compo- 
nent 1.2 magnitudes fainter at a distance of 1.79" (Fabricius 
& Makarov 2000). Thus, photoelectric photometric observa- 
tions of the system will include the secondary, which con- 
tributes 1/4 of the total light. (This also means that the 
amplitudes of the photometric variations we discussed here 
must be increased by one third to give intrinsic values.) 
The Tycho-2 photometry of the two components (Fabri- 
cius & Makarov 2000) transformed to the standard John- 
son system by using the relations by Bessell (2000) gives: 
V = 8.72, (B — V) = 0.20 for HD 99563A and V = 
9.91, (B — V) = 0.285 for HD 99563B. These values suggest 
that the two components could be physically associated. 

The HIPPARCOS parallax of the system (ESA 1997) 
is 4.27 ± 2.02 mas. The galactic reddening law by Chen et 
al. (1998) then results in A v = 0.010 ± 0.002. Therefore one 

ii for 



obtains M v = l.9t° 8 4 for HD 99563A and M v = 3.1; 
HD 99563B. Because of the large uncertainty of the HIP- 
PARCOS parallax we can therefore not obtain a radius es- 
timate of HD 99563A with sufficient accuracy, but we can, 
using also the (B — V) estimated above, point out that HD 
99563B is located inside the 5 Scuti instability strip and 
may therefore also show pulsations. However, our check for 
low- frequency variability did not detect these. 

To estimate the effective temperature of HD 99563 we 
can apply the calibration by Moon & Dworetsky (1985), 
who use the Stromgren Hp index as a temperature indica- 
tor. None of the other Stromgren indices is suitable for basic 
parameter determination of chemically peculiar A stars be- 



cause of heavy line blanketing. The measured Hp value for 
HD 99563AB is 2.830 (Olsen & Perry 1984). With the V 
magnitude and (B — V) colour differences given above, com- 
bined with the standard relation by Crawford (1979), we can 
determine Hp = 2.844 for HD 99563 A and Hp = 2.785 for 
HD 99563B. The calibration by Moon & Dworetsky (1985) 
then gives T aS ~ 8050 K for HD 99563A and T cff ~ 7400 K 
for HD 99563B. 

Elkin et al. (2005) performed a spectral analysis of 
HD 99563, which resulted in T cS = 7700 K, log g = 4.2 
and [M/H] = 0.5 for the roAp star. We assume a generous 
T of f = 7900 + 300 K and that HD 99563A is still on the main 
sequence, i.e. log g — 4.15 ± 0.2. Consequently, we estimate 
that R = 1.9l~o4 Rq. The projected rotational velocity of 
HD 99563A is 28.5 ± 1.1 kms" 1 (Elkin et al. 2005). Com- 
bined with the rotation period determined before, we then 
derive i = 601"^ degrees, a poor constraint. However, this 
also implies R ^ 1.58 Hq. 

As a final attempt to derive the basic parameters of HD 
99563A, we assumed that HD 99563B is a physical compan- 
ion. Due to the wide separation of the two components, they 
would have evolved independently from each other. The goal 
now is to find two stellar models that have effective tem- 
peratures as determined above (for reasons of consistency, 
we take the photometric values for both stars), that have 
the same age, and the observed magnitude difference of HD 
99563A and B. 

We used the Warsaw-New Jersey stellar evolution code 
(e.g. see Pamyatnykh et al. 1998 for a description) to find 
such a pair of models. We computed stellar evolutionary 
tracks in the range between 1.5 to 2.2 Mq for Z = 0.02. 
We indeed found two models that satisfy the observational 
constraints. They have masses of 2.03 Mq (HD 99563A) 
and 1.585 Mq (HD 99563B), luminosities of 21.6 and 7.1 L© 
respectively, and ages of 620 Myr. The model for HD 99563A 
has a radius of 2.38 Rq, which then implies i = 43.6 ± 2.1°. 
The internal accuracy of this method is high: changing the 
secondary mass by only iO.OlM© already cannot reproduce 
the observations within the errors. 

For the purpose of further discussion, we assume i — 
44°. We can now apply the OPM, which predicts 



An 



tan i tan f3 



(1) 



(Shibahashi 1986), where the A r 's are the amplitudes of 
the r th rotational sidelobes, i is the inclination of the stel- 
lar rotational axis to the line of sight and f3 is the mag- 
netic obliquity. From the amplitudes in Table 4, we obtain 
tan i tan/3 = 15.5 ± 1.5, hence (3 = 86.4 ± 0.3°. This is a 
perfectly reasonable result given the relative amplitudes in 
the central frequency triplet. 

Within the framework of a regular perturbation treat- 
ment, Shibahashi & Takata (1993) derived an expected am- 
plitude ratio 



^4-3 + A 3 _ 2. tan i tan /3, 
A2 + A-2 6 



(2) 



which is consistent with our observational data that lead to 
a null results when searching for the V2 and V-2 components 
of the stellar pulsation mode. 

We further examine the distorted dipole mode of HD 
99563 by applying the axisymmetric spherical harmonic de- 
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Table 5. Components of the spherical harmonic series description 
of the pulsation mode of HD 99563 for i = 44°, = 86.4°. 



1 





1 


2 


3 


(£) 

A_g [mmag] 








0.450 


A_ 2 [mmag] 






0.158 


0.158 


(£) 

A_[ [mmag] 




3.405 


0.037 


-0.699 


A\ [mmag] 


0.237 


0.376 


-0.093 


0.042 


A+i [mmag] 




2.372 


0.028 


-0.524 


^-+2 [ mma g] 
[mmag] 






0.088 


0.088 
0.180 


4>W [rad] 


0.024 


2.859 


-0.322 


2.820 


C ne n/K ma 3 




0.107 


0.020 


0.010 



Table 6. Comparison of the observed amplitudes and phases of 
the distorted dipole pulsation mode of HD 99563 and our spherical 



harmonic fit. 


ID 


A bs 




4>obs 


4>calc 


V 


0.29 ±0.03 


0.29 


+2.63 ± 0.09 


+2.63 




2.67 ±0.03 


2.67 


+2.87 ±0.01 


+2.87 


v +l 


1.82 ±0.03 


1.82 


+2.83 ± 0.01 


+2.87 


V-1 


0.00 (assumed) 


0.00 




n/a 


V+2 


0.00 (assumed) 


0.00 




n/a 


V-Z 


0.45 ±0.03 


0.45 


+2.82 ± 0.06 


+2.82 


V+3 


0.18 ±0.03 


0.18 


+3.00 ± 0.15 


+2.82 



composition method by Kurtz (1992) which is based on the 
theory by Shibahashi & Takata (1993), to our data. This 
technique breaks the magnetically distorted mode up into 
its pure I — 0,1,2,... spherical harmonic components and 
consequently allows one to infer the shape of the mode. We 
applied the method to our frequency solution (Table 4), as- 
suming zero amplitude for the unobserved compo- 
nents of the mode. The result is given in Tables 5 and 6. 

From the values in Table 5 it can be seen that the mode 
of HD 99563 is dominated by the dipole component, which 
is reasonable as the magnetic fields in Ap stars are predom- 
inantly dipoles. The contribution of the £ — and £ = 2 
terms to the mode are fairly small, but the I = 3 com- 
ponent does have some influence. Table 6 shows that the 
observations are quite well reproduced by our fit. It would 
be interesting to see whether the theory by Saio & Gautschy 
(2004) is capable of reproducing our results. 

We can now check how well our model reproduces the 
observed pulsational amplitudes and phases over the stellar 
rotation cycle. To this end, we subdivided the time series 
into pieces some 4 — 5 pulsation cycles long and determined 
the amplitudes and phases for the 1557.653 /iHz variation 
within these subsets. The spherical harmonic decomposition 
method by Kurtz (1992) yields a fit to the amplitude/phase 
behaviour over the rotation cycle, and we show it together 
with the data in Fig. 6. 

We can see that both magnetic (pulsation) poles are 
actually in view for approximately the same amount of time; 
the relative fractions are 53 and 47 per cent. This is also 
the reason why it was so difficult to determine the rotation 




Rotation Phase 



Figure 6. Pulsational amplitudes and phases relative to a rota- 
tion period of 2.91179 d. Phase zero corresponds to pulsation am- 
plitude maximum, and one and a half rotation cycles arc shown. 
The line is a fit computed with Kurtz' (1992) spherical harmonic 
decomposition method. 

period of the star; the pulsations as well as the mean light 
variations are almost symmetrical along the rotation period. 

The fitted curves reproduce the observations fairly well, 
with one exception: the pulsation amplitude is overesti- 
mated for the pole that is for a shorter time in the line 
of sight, whereas the amplitude is underestimated for the 
other pole. The harmonic frequencies cannot be made re- 
sponsible for this, because they do not affect the first-order 
amplitudes analysed here. Consequently, we believe that the 
poor fits near the pulsation amplitude maxima could be due 
to presently undetected additional multiplet components of 
the pulsation mode; the amount of over- and under-fitting 
can be fully explained with the noise level in our residual 
amplitude spectrum. 



7 CONCLUSIONS 

Our photometric multisite observations of the roAp star HD 
99563 resulted in the detection of a frequency quintuplet 
that is due to a single distorted dipole pulsation mode. The 
splitting within this frequency quintuplet together with our 
mean light observations and published magnetic measure- 
ments allowed us to determine the stellar rotation period as 
2.91179 ± 0.00007 d. Within the errors, the mean light ex- 
trema occur in phase with pulsation amplitude maximum, 
suggesting that the abundance spots on HD 99563 should 
be fairly concentric around its magnetic poles. 

To our knowledge, HD 99563 is only the fourth roAp 
star where both magnetic poles become visible throughout 
a rotation cycle (the others are HR 3831, e.g., see Kurtz et 
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al. (1993), HD 6532, e.g., Kurtz et al. (1996), and HD 80316 
e.g., Kurtz et al. (1997)), and only the third whose pulsa- 
tional mode distortion has been quantified. In this context it 
is interesting to note that the observational features of HD 
99563 are in many aspects similar to another roAp star, the 
well-studied HR 3831 (e.g., Kurtz et al. 1993, Kochukhov 
2004): their geometrical orientations, rotation period and 
effective temperatures are alike. Only the pulsation period 
of HD 99563 is some 10 per cent shorter, and the phases its 
first-order combination frequencies are more consistent with 
them being harmonics. It would thus be very interesting to 
compare these two objects in detail. 

This would however require more in-depth studies of 
HD 99563. Its basic parameters are still poorly known; the 
radius we inferred is based on the assumption that the visual 
companion is physical, which needs to be checked. If HD 
99563 A and B were a physical pair, then the separation of 
the components is too wide for motion around a common 
centre of mass to be determined in a reasonable period of 
time. However, common proper motion may be detectable 
within a few years. 

The most efficient way to pin down the radius of the 
roAp star would perhaps come from a combination of mag- 
netic and polarisation measurements (Landolfi et al. 1997). 
With that, the inclination i and the magnetic obliquity f3 
can be constrained, and given the accurate rotation period 
and v sin i already available, a fairly accurate radius could 
be obtained. 

There are several other possibilities that make HD 
99563 interesting for future observations. As already ar- 
gued by Elkin et al. (2005), the star is a very attractive 
target for pulsational radial velocity measurements. Given 
the high amplitude of the mean light variations of the star 
and its favourable geometry, Doppler Imaging of its surface 
should also be within reach. Finally, we have evidence that 
we have not yet deciphered the full pulsational content of 
the star's light curves. Another multi-site campaign, aim- 
ing at obtaining > 200 h of observations on 1-metre-class 
telescopes would therefore also be justifiable. As we pointed 
out, HD 99563B is located within the 5 Scuti instability strip 
and should therefore also be tested for pulsations, with more 
suitable means than ours. 
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